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Abstract 

We propose a non-thermal scenario for the generation of baryon number asymmetry in a 
radiative neutrino mass model which is modihed to realize inflation at the early Universe. 
In this scenario, inflaton plays a crucial role in both generation of neutrino masses and 
lepton number asymmetry. Lepton number asymmetry is hrstly generated in the dark 
matter sector through direct decay of inflaton. It is transferred to the lepton sector via 
the dark matter annihilation and then converted to the baryon number asymmetry due 
to the sphaleron interaction. All of the neutrino masses, the baryon number asymmetry 
and the dark matter are intimately connected to each other through the inflaton. 
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1 Introduction 


Recent experimental and observational data for neutrino masses DE] and dark matter 
(DM) |3l m |5] suggest that the standard model (SM) should be extended. The radiative 
neutrino mass model proposed in [B] is such a simple extension of the SM with an inert 
doublet scalar and right-handed neutrinos. It seems to be a promising candidate which 
could take the place of the famous canonical seesaw model for neutrino masses [7]. An 
interesting point of this model is that it could also give the origin of DM [HI |9]. A Z 2 
symmetry imposed to forbid the neutrino masses at tree-level could guarantee the stability 
of the lightest Z 2 odd held, which could be DM. In this model, DM is an indispensable 
ingredient for the neutrino mass generation at TeV regions. 

Although the model has such interesting aspects, baryon number asymmetry in the 
Universe [10], which is another crucial problem of the SM, cannot be easily explained in 
a consistent way with the relic abundance of DM. If we suppose the ordinary thermal 
leptogenesis mm, the sufficient baryon number asymmetry can be generated only in 
the case where the model has a hnely tuned spectrum for the Z 2 odd helds. 

If the lightest right-handed neutrino is assumed to be DM, both its relic abundance and 
small neutrino masses require 0(1) neutrino Yukawa couplings in generajl] [8|. They can 
allow to cause large CP asymmetry in the decay of right-handed neutrinos even if their 
masses are of 0(1) TeV. However, the same neutrino Yukawa couplings could cause large 
washout of the generated lepton number asymmetry through the inverse decay and the 
lepton number violating scattering processes. As a result, the thermal leptogenesis is not 
easy to generate sufficient lepton number asymmetry in a consistent way with the neutrino 
oscillation data and the DM abundance at least in the simplest form of the model [13] . On 
the other hand, if the lightest neutral component of the inert doublet scalar is assumed to 
be DM [13], the neutrino Yukawa couplings could be small enough to be consistent with 
both the DM relic abundance and the small neutrino masses. However, the large CP 
asymmetry in the decay of right-handed neutrinos requires fine mass degeneracy among 
the right-handed neutrinos [I5]. Non-thermal leptogenesis m [H] might give another 
consistent scenario for the origin of the baryon number asymmetry in this model or its 
supersymmetric extension |18] . 

^ This brings about dangerous lepton number violating processes at large rate unless special flavor 
structure is assumed for the neutrino Yukawa couplings [Sj. 
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In this paper, to solve the above mentioned fault for leptogenesis, we propose a simple 
scenario in the model which is extended so as to incorporate the inflation at the early 
Universe [19]. The neutrino mass generation is connected with the inflation through the 
inflaton interaction. The lepton number asymmetry is also produced through the inflaton 
decay in the inert doublet sector which contains the DM candidate mi OS]. After this 
lepton number asymmetry is transferred to the lepton sector via lepton number conserving 
scattering processes, the sphaleron interaction converts a part of it to the baryon number 
asymmetry. 

Remaining parts of the paper are organized as follows. In the next section, we intro¬ 
duce the extended model briefly. In section 3, we study its phenomenological features. 
Firstly, we describe the inflation in the model and also the small neutrino mass gen¬ 
eration. After that, we explain the scenario for the generation of the lepton number 
asymmetry and then estimate the baryon number asymmetry expected to be produced 
hnally. Following this discussion, the consistency of the scenario with DM phenomenology 
is examined. Relation between the present DM scenario and the asymmetric DM scenario 
is also remarked. We summarize the paper in section 4. 

2 An extension of the radiative seesaw model 

Our model considered here is based on the one proposed for the radiative neutrino mass 
generation [6]. The original model is a simple extension of the SM with an inert doublet 
scalar r] and three right-handed neutrinos N^.. These new helds are assigned odd parity 
of an imposed Z 2 symmetry, although all the SM contents are assumed to have its even 
parity. Invariant Yukawa couplings and scalar potential which are relevant to these new 
helds are summarized as 

- = h.,NK,ri'iL. + htf(L.i)NR, + i [MtNR.N'k, + MiN^NR,) , 

-I- + rri^V^T] + Ai(0V)^ + ^2i'n^vy + 

+ y + h.c.] , (1) 

where iii is a left-handed doublet lepton and (j) is an ordinary doublet Higgs scalar. We 
use the basis for which both matrices for charged lepton Yukawa couplings and right- 
handed neutrino masses are real and diagonal. Since the Z 2 is assumed to be the exact 
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symmetry of the model, the new doublet scalar rj should not have a vacuum expectation 
value. As its result, neutrino masses are forbidden at tree level and the lightest held with 
the odd parity is stable to be DM. 

In this type of model, the lepton number L is usually assigned to these new helds as 
Lij]) = 0 and L{Nr.) = 1. In such a case, the neutrino mass generation and leptogenesis 
have been studied under the assumption that mass terms of the right-handed neutrinos 
violate the lepton number [131 [IS] • The DM abundance has also been studied supposing 
that either the lightest right-handed neutrino or the lightest neutral component of rj is 
DM. However, it is useful to note that there could be another assignment of the lepton 
number such as L{r]) = 1 and L{Nr.) = 0 [1^. In this case, is forbidden as long 

as the lepton number is imposed as the exact symmetry. As a result, neutrino masses 
could not be generated even if the radiative effect is taken into account. Thus, some 
suitable origin of the lepton number violation should bring about this As term as an 
effective interaction at low energy regions. We study such a possibility in the following 
part. 

For this purpose, we consider an extension of the model at high energy regions by 
introducing canonically normalized complex singlet scalars Sa which are assigned odd 
parity of the Z 2 symmetry and L = 1. The potential and interaction terms of are 
assumed to be given by 
2 

~ Cs = + K3{SlSa){r]^r]) 

a=l 


+ 

+ 



.SlSi 


Is? 




( 2 ) 


where both n and m in the third line are positive integers and Mpi is the reduced Planck 
mass. Although the Z 2 is kept as the symmetry of these terms, the lepton number is 
violated through the mass terms in the second line and also the Planck 

suppressed C 2 terms in the third line. The latter one is neglected in the low energy region. 
On the other hand, the former lepton number violation could be an origin of A 5 term in 
eq. (II]). In fact, as a simplest case, we might consider the situation where ^ is 


4 









satisfied. In this case, the model defined by eq. (II]) can be easily obtained as the effective 
one with A 5 = Hq, where Ag^^ is defined by Ag^^ = They are induced as the 

Sot 

effective interaction terms at low energy regions after the singlet scalars are integrated 

out [T71 [19] . 

In the following discussion, we are focus our study on the situation such that the terms 
in the last line in eq. (E]) could be a dominant part of the potential at the early Universe. 
We suppose that |S'i| takes a large but sub-Planckian value in such a period. It could be 
realized under the condition such aj^ 


Ki C Cl 





2 


5 




2 


< Cl 



(3) 


where ipi is defined by Si = ^ 6 *^^ and ipi < M^\. If we use the polar coordinate of Si 
defined here, the last line of eq. ([ 2 ]) can be written as 


= Cl 




2nJp^2n-4 


1 + 2C2 




2m 


COS 


2 A 2 


+ 2m9i 


,a/2Mpi 

We easily hnd that has local minima with the potential barrier 14 ~ 


(4) 


2(n + 7Ti) 


2^1 + 71 




in the radial direction, which form a spiral-like trajectory. We consider the inflation which 
is caused by the inflaton evolution along this trajectory. 


3 Phenomenological features of the model 

3.1 Inflation 


We briefly review the features of the inflation induced by the potential 01]). We assume 
that ifi takes a large initial value on a local minimum in the radial direction. In that 
case, as shown in [19], the model could cause sufficient e-foldings through the inflaton 
evolution along the spiral-like trajectory even for sub-Planckian values of (pi. An inflaton 
field X could be identified with 

^ 6mA2 6mA2 ’ 

^When Si plays a role of inflaton, this condition could be relevant to the rj problem in this inflation 
scenario. We cannot fix it at this stage unless the UV completion of the model is clarified. 
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Cl 

C2 

A 

Mpi 

Fi 

\/2Vfpi 

H, 

N, 

ris 

r 

(xlO-^) 




(xlO^^GeV) 




9.84 

1.7 

0.05 

0.411 

5.91 

60.0 

0.964 

0.056 

8.62 

1.9 

0.05 

0.406 

5.40 

60.0 

0.959 

0.040 


Table 1. Examples of the predicted values for the spectral index Us and the tensor-to-scalar ratio r in 
this scenario fixed by n = 3 and m = 1. 


where the field a is defined as 


da = 


(pi + 




1 + 4m^ 



4 l 1/2 

ipidOi- 


( 6 ) 


Fields with the subscript e stand for the fields at the end of inflation. The number of 
e-foldings caused by x is given as 


N = 


1 


rxe 


V, 


dXTfr^N{x)-N{Xe), 


Jx 


Vl 


(7) 


Si 


where and N{x) is represented by using the hypergeometric function F as 

1 


N{x) = 


QnFn \ A ) V\/2Mpi 


^ 6c2m / (fi 

n{3 + m) V\/2Mpi 


2m 




2 m' 


( 8 ) 


Here we note that the model could have a different feature from the ordinary inflation 
scenario such as the chaotic inflation. In eq. ([7j), N{x) Af{xe) might not be satisfied 
generally. In this model, inflation is expected to end at the time when — Id, is satisfied. 
If we apply the slow-roll approximation 3Hx = ~'^Si slow-roll parameters 

e = inflation is found to end at e = This means that the end 

of inflation could happen much before the time when e ~ 1 is realized since > Vb 
is satisfied. In that case, N{xe) could have a substantial contribution to determine the 
e-foldings N in eq. O- 

The slow-roll parameters e and r] = Mpj (j can be represented by using the model 
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parameters as 


e = 




2 


rj = m 


f V^MpA " / ^ ^(2^ - 3) - 2c2(m + n)(2m + 2n - 3) (7^) 

I 71 ] VMpi 


2m 


1-2c2 


\ 2m 


\/2Mpi 




( 9 ) 


If C 2 terms are neglected in these formulas, we find very simple formulas for these slow-roll 
parameters at the period characterized by the inflaton value y*. They can be represented 
by using the e-foldings dehned for N{x*) in ^q. (IHl) as 

2n — 3 


n 


e ~ 


6 (iV, + iV(y,))’ 6{N, + N{xe)) 

Thus, the scalar spectral index n* and the tensor-to-scalar ratio r can be derived as [19 

n + 3 _ 8n 


( 10 ) 


Us = 1 — 6e + 2?/ ~ 1 


3(iV, + iV(x,))’ 


r = 16e ~ 


3(iV, + iV(ye)) 


( 11 ) 


If we focus on the case n = 3, these formulas reduce to the ones of the chaotic 

inflation scenario [21]. However, as shown in [19], the values of rig and r in this model 
could deviate from the ones of the chaotic inflation due to the non-negligible C 2 

term contribution. Taking account of uncertainty caused by the reheating process and 
others, might be considered to take a value in the range 50 - 60. If we estimate both 
rig and r by hxing the parameters in the potential suitably, they could take consistent 
values for in this range with the ones suggested by a joint analysis of BICEP2, Keck 
Array and Planck [22l [23] . Such examples for n = 3 are shown in Table 1. The condition 
(]3|) requires fhs^ -C 10^^ GeV in this case. Much better agreement with the observational 
results for Ug and r is found in the case n = 1,2 [19] . 

Finally, we note that the polar coordinate cannot be used for to rewrite the potential 
as eq. (|1|) unless = 0 is satished. In order to make this inflation scenario possible, 
should be generated after the end of inflation at least. It is not difficult to modify the 
model to satisfy this condition. For example, we may introduce a singlet scalar tf) with 
L = —1. In this case, its potential might be given by 


Vy = + {C 2 SIS 1 - + h.c.). (12) 
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If the value of |S'i| becomes smaller than y after the end of slow-roll inflation, -i/; could 
get the vacuum expectation value which induces the required mass term for Sa through 
the ^3 term. After the generation of these terms in eq. ([ 2 ]) as the effective ones, the mass 
splitting between the real and imaginary components of Sa is brought about. Each mass 
eigenvalue is expressed as ± where -|- and — signs correspond to the 

real and imaginary component, respectively. We note that the stability of the vacuum 
requires > nT,|^. The difference of these mass eigenvalues can be a measure of the 
lepton number violation in the model. 


3.2 Neutrino masses 


The neutrino masses are generated in the similar way to the original model. The one- 
loop effect which picks up the lepton number violation induced by the mass term 
generates the neutrino masses through the electroweak symmetry breaking as shown in 
the left-hand diagram of Fig. 1. The neutrino mass matrix obtained in this way can be 
described by the formula 

(•^ “ZEE A4. m,„). (13) 

k=l a=l,2 f=± 


where = 171 ^ + (A 3 -|- and (0) = 174 GeV. stands for ^ and 

respectively. The function I{ma,mb,mc) is dehned as 


I{ma,mb,mc) = 


+ 


(m^ — Inm^ ml Inm^ 

{ml — miy{ml. — mlY {ml — ml){ml — mlY 
ml In ml 1 

{ml — ml){ml — mlY {ml — ml){ml — ml) 


As long as is satished, this formula is found to be reduced to 


(14) 


M 


u 

st 


hskhtk{4>)‘^ 

k=i ^ 




(15) 


where we neglect logarithmic factors. If we note that two right-handed neutrinos are 
enough to explain the neutrino oscillation data, hi could be assumed to be so small that 
the contribution of Ni to the neutrino masses is negligible. We adopt this assumption 
throughout the following discussion, for simplicity. 
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Fig. 1 Left: a one-loop diagram contributing to the neutrino mass generation. The dimensionful 
coupling is defined as ^ and /ii ^ by using fia in eq. ([2|). Right: a one-loop diagram 

contributing to the lepton flavor violating process £i —> ijj. 


If we assume the flavor structure of the neutrino Yukawa couplings discussed in Ap¬ 
pendix A, the required mass difference for the atmospheric neutrinos and the solar neu¬ 
trinos could be explained by the largest mass eigenvalue and the next one in this mass 
matrix, respectiveljo. For example, this requirement could be represented as 

1-2 \ 2 


Y 

a=l,2 

Y 

0 = 1,2 


he 


ho 


mz 


ho 


m: 


- 1-0 


ho 


~ 10 


-6 


m: 


+a. 


~ 10 


-6 


5.1 X 10- 

h2 

2.7 X 10-2 
h.s 


Mo 




2 X lO^GeV J ’ 

M 3 \ 

5 X lO^GeV ) ’ 


(16) 


where we assume Mn = 1 TeV and CP phases are neglected in this estimation. It should 
be noted that the left-hand side of eq. ffTOl) corresponds to the effective coupling A 5 . It 
plays a crucial role also in the generation of baryon number asymmetry and DM direct 
search as discussed later. 

It is well-known that these new helds induce the lepton flavor violating processes at 
one-loop level. The typical one is £i —>■ £j^ whose diagram is shown in the right-hand side 
of Fig. 1. Its branching ratio can be estimated as m 




So; 


647r(Gi.M2)2 


^ ^ hikhjkF2 




k=l 


Vm, 


~ 8 X 10-^ 


hikhjkF2 ( ^ 


k=l 


(17) 


^It should be noted that one of the eigenvalues of this assumed mass matrix is zero. It may be also 
useful to recall that the cosmological upper bound for the neutrino masses is 0.23 eV [2^ . 
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( 18 ) 


where 


1 TeV is used and F 2 {x) is given by 


F2{x) 


1 — + 2 x® — In 

6(1 — x‘^Y 


Here we note that F 2 {x) ~ ^ for x 1 and the present upper bounds for Br{fi —)■ ey) 
and Br{T fi'j) are given as 5.7 x 10“^^ [25] and 4.4 x 10“® [26], respectively. Since 
Mfc > Mri is assumed in the present model, the bounds for these flavor violating processes 
give no substantial constraint on neutrino Yukawa couplings as found from eqs. ffT6D and 

(ED. 


3.3 Baryon number asymmetry 


Reheating process should follow the inflation discussed in the previous section. In this 
scenario, reheating is expected to occur through the decay of Si after the inflaton stops 
its evolution along the above mentioned spiral-like trajectory and S'±i starts to oscillate 
around a global minimum of the potential. Although preheating could occur via scalar 
quartic couplings in the hrst line of eq. ([ 2 D, the reheating is expected to be hnally com¬ 
pleted through the decay of Si [27ll2H|- Since lepton number asymmetry is not produced 
through the particle creation in the preheating, we focus our study on the decay of Si 
here. 

The decay of Si is induced by the interaction of Si with 0 and rj during the oscillation 
induced by the mass terms which are given in the second line of eq. ([2D. The reheating 
temperature may be estimated by using the usual instantaneous thermalization approxi¬ 
mation. If we use this approximation, the reheating temperature is determined through 
the condition if ~ T^i. if is the Hubble parameter and r±i stands for the decay width of 
Ail which is the real and imaginary component of Si. Since r±i can be approximately es¬ 
timated as r±i ~ where = ifig^ ±m|^, the decay products of S±i —)■ r](j)\ 

hnally make thermal plasma with possible reheating temperature 


ry>~0.35jX'"‘lftl (^)’. (19) 

\m±ij 

where we use = 116 as the relativistic degrees of freedom in this model. If we consider 
a situation such that S±a is not thermally generated through the inverse decay or the 

^ In this estimation, the oscillation energy of each component is assumed to dominate the total energy 
density of the Universe. 
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scatterings, m±a > should be satisfied at least. This condition could be expressed as 


/ii 

m+i 


< 1.9 X 10 


-4 


fe)( 


m+i \ 2 


\m+J U09 GeV/ 


( 20 ) 


In the following part, we confine our study to the case where this condition is satisfied. 

The infiaton decay is relevant to the generation of baryon number asymmetry in this 
model. The lepton number asymmetry could be directly generated through this process 
non-thermally since this decay violates the lepton number. In fact, if is complex, the 
cross term between tree and one-loop diagrams for the decay could bring about the CP 
asymmetry. The CP asymmetry induced through this decay of S'±i can be estimated aj^ 


e± = 


T{S±i -)■ 7701) -h T{S±i -)■ r]^) 


|/72p sin 2(01 - 02 ) 


IbTT 


(m|i + 

m^i {jn\i + ?77^2)^+2 


In 
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±1 


m 


+2 


— m 


2 

-2 


[m 


±1 


m 


+2 


)2 + ( 


m 


±1 


m 




( 21 ) 


where 0* = arg(pj) and r±Q, = ^1 — -^r-^ ■ As long as the condition (1201) is satisfied, 

the lepton number asymmetry generated through the infiaton decay could be the only 
source for the baryon number asymmetry since there is no mother particles in the 
thermal bath. 

If both components S±i have finely degenerate masses ~ their decay occurs 
almost simultaneously and then ~ We also find that e+ ~ —e_ is satisfied. 

Since the lepton number asymmetry generated in the rj sector through this decay could 
be estimated as AL ~ e+n 5 ^^(T^’^^) + e-ns_i{T^ ^), AL may not take a large value in this 
case because of the cancellation due to e-ns_^ (T^ ^) — —e+nsj^^ On the other hand, 

if substantial mass splitting appears between the components and then is 

satisfied, the S'+i decay is expected to occur later compared with the decay of S_i because 
of r_i > T+i. In such a case, a part of lepton number asymmetry generated by the S-i 
decay could be washed out by the lepton number violating processes before the delayed 
S'+i decay. Thus, the lepton number asymmetry expected in the r] sector after the S'+i 
decay could be estimated as AL ~ e+n 5 _^^(T^^^) -|-/C^(T^^^)e_n 5 _^(r^ ^) where /Ci„(T^^^) 
represents the washout effects from ^ to If the lepton number violating processes 


^In the following study, we assume the maximum CP phase | sin 2 ( 0 i — 6 * 2 )| = 1 - 
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Fig. 2 Feynman diagrams which contribute to the transfer and the washout of the lepton number 
asymmetry. The left diagrams are lepton number conserving scattering processes whose reaction densities 
are represented by 70 (upper ones) and (lower one). The right diagrams are lepton number violating 
scattering processes whose reaction densities are represented by (upper ones) and 7 j, (lower one), 
respectively. 


decouple and then /C^ = 1 is satisfied in this period, AL is expected to take a substantial 
value because e_ns'_^(T^~^) 7 ^ —(T^"^^) is satisfied. 

The lepton number asymmetry generated in the i] sector via the S±i decay cannot 
be transferred to the SM contents through the decay of rj. We should note that 7 ] does 
not have any decay modes to the SM contents because of the Z 2 symmetry. However, it 
could be partially transferred to the lepton sector through the lepton number conserving 
scatterings 7777 ii and rji —)■ rjH. These are induced by neutrino Yukawa couplings 
and their diagrams are given in the left-hand side of Fig. 2. On the other hand, it could 
also be washed out through the lepton number violating scattering processes 7777 —)■ ( fxj ) 
and 7701 —)■ rf'cj). These are caused by the S±a exchange due to the 77 ^ couplings. Their 
diagrams are also shown in the right-hand side of Fig. 2. In the situation where these 
processes are competing with each other before reaching the weak scale, the lepton number 
asymmetry kept in the lepton sector could be converted to the baryon number asymmetry 
through the sphaleron interaction. We examine this scenario quantitatively by solving 
relevant Boltzmann equations. 

For this purpose, we define the lepton number asymmetry in the co-moving volume 
as AW = AAh the lepton sector and AW? = s ^ sector, respectively. 

The entropy density s is expressed as s = As discussed in the previous part, the 

lepton number asymmetry in the 77 sector is expected to be fixed through the decay of S'±i. 
Thus, at the reheating temperature the lepton number asymmetry in each sector 

are supposed to be AY£(T^’^^) = 0 and AW?(T^^^) = ^ where sr 
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stands for the entropy density at If we use ^ and Ps±i{Tj^^) 

which are derived by assuming the instantaneous thermalization after the S. 
decay, we hnd that the latter can be expressed as 


±1 


q rp{+) 


3 T, 


(-) 


+ -e--^ 
4 m_i 


( 22 ) 


By taking account of the relevant processes which are explained above, Boltzmann 
equations which describe the evolution of All and Al^ are given a^ 


dAYr, 

dz 

dAYe 


sH{M^ 


2(7a + 7b) 


AyA , ,AK 

+ 2(73; + 7 y ) 


Yr 




-2(7a + 76) 


/AA Ay 


Ur 


yr 


(23) 


in the 


dz sH{My) 

Since we consider the case where the condition (l 2 n|l is satished, the effect of S±a 
thermal bath can be neglected. Each reaction density 7 * is explained in the caption 
of Fig. 2 and their formulas are given in Appendix B. The generated baryon number 
asymmetry could be estimated as irzi 


Yb — ——^Yi{zEw) 


(24) 


by using the lepton number asymmetry Ay obtained as the solution of these equations 
at the weak scale. 

Although detailed analysis of the generated baryon number asymmetry requires to 
solve the above Boltzmann equations numerically, we briefly discuss their qualitative 
aspects before proceeding to it. At hrst, we note the behavior of the ratio of the re¬ 
action rate T to Hubble parameter H for the relevant scattering processes in the case 
m±a > which we consider here. T and H are expressed as 1^ 6 = 
where ~ — -^z~^K 2 {z) and H{z) ~ 0.33gy‘^^^z~‘^. In the lepton 

number conserving scattering processes caused by the neutrino Yukawa couplings, 
is a convex function of 2 ; which takes a maximum value around Zm — They freeze out 
at Zf{> Zm) in the case where > 1 is satished at Zm- It is important to note that 

Ay follows Ay to be Ay = Ay^ as long as ^ 1 is satished. On the other hand, 

the coupling which causes the lepton number violating scatterings is dimensionful so 

^Following the usual convention, we introduce a dimensionless parameter z as z = ^ by using a 
convenient mass scale My, which is defined below eq. (HI. 
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r H-F M 

that ""jj ^ increases monotonically with z throughout the range —^ < z < 1. Since 

R 

these processes are expected to be in the thermal equilibrium at a certain period Zg where 
= 1 is satisfied, is expected to be erased at z ^ Zg. However, these processes 
are suppressed at z ^ 1 by the Boltzmann factor. 

Here we note that both Zf and Zg are determined by the parameters relevant to the 
neutrino masses. We could make a rough estimation of favored parameters for the gen¬ 
eration of baryon number asymmetry by taking account of it and the above arguments. 
As seen in eq. ffTbjl . the neutrino oscillation data imposes a relation for neutrino Yukawa 
couplings and a GeV unit such that 


{hh^)kk 

Mk 
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-ho: 


0(10 




(25) 


where we assume = 1 TeV. If we use this condition, both Zf and Zg can be roughly 
estimated as 


z, ~ O(10‘«) 


Ml 


0(10 


-13n 


. 0 = 1,2 




m_ 


0(10 


ho 


-111 
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+a 




ho 


m_ 


-2 



(26) 


where the CP phases of neutrino Yukawa couplings are neglected. These results suggest 
that Zf > Zg is always satisfied. 

The washout factor /Cu,(z) which we have already introduced in the previous discussion 
is characterized as a decreasing function at z > Zg and /C^(z) ~ 1 at z < If we use it, 
the total lepton number at z might be written as 


AW(z) + AY,(z) = /C^(z)AY,(zk), (27) 

where we use eq. fl22|) as the initially generated lepton number asymmetry. On the other 
hand, the lepton number asymmetry in both sector at z could be related as 


AW(z) =/Ct(z)AY^(z), (28) 

where /Ct(z) stands for the transfer efficiency of the lepton number asymmetry from the r] 
sector to the doublet lepton sector. If the lepton number conserving scattering processes 
are in the thermal equilibrium, /Ct(z) = 1 is satisfied. Using these relations, we could 
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h2 


hs 

msi 

mSi 

Imi| 

^S2 

ms2 

ms. 

IM2| 

ms. 


e+ 


(a) 

1.0 

■ 10"2 

4.8 

■ 10"3 

0.5 

2.0 • 10"® 

1.3 

0.5 

3.0 • 10-3 

1.7 

• 10-3 

1.0 • 10-^° 

(b) 

1.8 

• 10"2 

9.5 

• 10"3 

0.5 

10"® 

1.3 

0.5 

8.0 • 10-2 

1.2 

• 10-2 

3.1 • 10-® 


Table 2. The CP asymmetry e+ and the baryon number asymmetry \ Yb \ obtained in the present scenario 
for typical parameter settings. The dimensionful model parameters are taken to be (a) M 2 = 2 x 10“^, 
M3 = 5 X 10^ and ifiSi = 10®, (b) M2 = 2 x 10®, M3 = 5 x 10® and rhsj = 10® in a GeV unit, respectively. 
Neutrino Yukawa couplings are numerically determined for M^ = 1 TeV so as to realize the neutrino 
mass eigenvalues required from the neutrino oscillation data. 

consider two possible cases for the generation of lepton nnmber asymmetry in the lepton 
sector. 

(a) If the lepton nnmber conserving scatterings are in the thermal eqnilibrium at an 
early stage and freeze ont at Zf, the lepton number asymmetry in the lepton sector at the 
weak scale is found to be roughly expressed as 

^ Aiy.-H). (29) 

1 + ICt{zf) 

Although lCw{zf) = 1 is satished for Zf < z^, the neutrino mass condition allows only the 
situation Zf > Zf, as shown in eq. fl26l) . Thus, the required value of ^Y^^zew) could be 
obtained in the case where K^{zf) is not so small. It could be realized only for S> M^j. 

(b) If the lepton number conserving scattering processes never reach the thermal equi¬ 
librium at z{< Zg) but has non-negligible values, the situation becomes completely 

different from the case (a). In this case, a part of could be transferred to the lep¬ 
ton sector. Since steeply decreases at z ~ Zg, AYg could take a hxed value which 
might be roughly estimated as AYg(zg) independently of the value of at z(> Zg). 

The transferred lepton number asymmetry AYglzg) is kept until the weak scale. Thus, 
AYg[z-Ew) could be expressed as 

AI^(^ew) — ^t{.Zg) AY^{zr) . (30) 

where K,t{zg) <C 1. Thus, the required lepton number asymmetry in the lepton sector 
could be obtained at the weak scale for a suitable Kt{zg). Such a situation could happen 
only in the case Mk ^ 

Now we present results of the numerical analysis of the Boltzmann equations. Model 
parameters used in this analysis are summarized in Table 2, which are numerically fixed 
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Fig. 3 The left-hand panels show the results for the case (a). The ratio of reaction rate T to the Hubble 
parameter H for each relevant process is plotted as functions of z in the upper panel. The solutions 
AYri and AYg of the Boltzmann equations are shown as functions of z in the lower panel. The lepton 
number asymmetry required to explain the observational results is shown by the horizontal black line. 
The right-hand panels show the results for the case (b) in the same way as the case (a). 

to satisfy the conditions for the nentrino masses. If we take acconnt of the conditions 
ffT^ and fl 20 l) . we find that shonld be satished and also their phases can be 

hxed as 01 7 ^ 0 and 02 = 0, |. This justifies the estimation in eq. ffTbjl . fl25l) and fl26|) and 
the assumption for the maximum CP phase in eq. (El]), which is used in this analysis. It 
also allows r]R and rji to be the mass eigenstates of the neutral components of rj. This 
becomes important for the study of DM phenomenology in the next subsection. 

Solutions of the Boltzmann equations fE3|) for these parameter settings are presented 
in Fig. 3. In the upper panels of this figure, the for the relevant processes are plotted as 
functions of z. In the lower panels, and AY^ are plotted as functions of z. The lepton 
number asymmetry required for the suitable baryon number asymmetry is also shown by 
the horizontal black dotted lines in these panels. The left and right panels show the 
results corresponding to the cases (a) and (b) discussed above, respectively. They show 
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that the above discussion describes qualitatively the features of the present scenario well. 
Although our study here is done only for the limited parameter sets, the results show that 
the scenario could generate the sufficient baryon number asymmetry for suitable model 
parameters in each case. Detailed study of this scenario for wider range of the model 
parameters will be given elsewhere. 

We should recall again that the same parameters used here are closely related to 
several low energy phenomena. Although some of them have been discussed already, 
there is another one which has not been taken into account still now. We need to check 
the consistency with it to see whether the model works well or not. It is DM physics and 
this issue is the subject in the next part. 

3.4 Dark matter 

The DM candidate is built in the model as the lightest Z 2 odd held. We identify it as the 
lightest neutral component of rj. We choose to be real and -C is supposed 
to be satished. In this case, the real and imaginary parts of the neutral component of r], 
which are written as r]ji and rji, become the mass eigenstates as mentioned before. If r]ji 
is supposed to be a DM candidate, r]R could be scattered with nuclei inelastically to rji. 
It is mediated by the Z boson exchange. Since it contributes to the DM direct search 
experiment [29], a strong constraint is imposed on the mass difference 5(= 
between r]R and ng. This might give the scenario an interesting chance for giving a 
prediction in the DM direct search experiments as seen below. 

We recall the experimental situation that we have no evidence in the DM direct search 
experiments [30]. If we apply it to the above mentioned process, we could put a bound 
for 5. It might be estimated as 5 > 150 keV conservatively. Since this mass difference is 
expressed in the present model as 



(31) 


the constraint is found to be represented as 



(32) 


"^The mass of rjn and rji can be expressed as and respectively, 

by using the effective coupling A 5 . 
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As noted in the previous part, the left-hand side of eq. (l32il corresponds to the effective 
coupling I As I for the assumed parameters. Although this constraint depends on the DM 
velocity distribution in our galaxy and other uncertain factors, eq. fl3^ gives an interesting 
condition for the present scenario on the origin of the baryon number asymmetry. We 
hnd that the model parameters used in the case (b) gives IA 5 I ~ 3 x 10“^ and then 
this condition is clearly satished. On the other hand, the situation is subtle in the case 
(a) since we hnd IA 5 I ~ 5 x 10“®. This suggests that the DM candidate in this model 
could be detected through the inelastic scattering in the direct search experiments if this 
leptogenesis scenario is realized in Nature for this parameter range. It may be worthy to 
reexamine the direct search results in this mass range in detail. 

The above scenario should be also consistent with the DM relic abundance. In the 
present study, DM is assumed to be rjR. In general, its relics could come from two types 
of origin such as 

VLh? = -|- Dnonth^^- (33) 

The hrst one is the usual thermal relic, that is, the remnant of rj^ decoupled from the 
thermal equilibrium distribution. It can be estimated by using the usual formulas |31] . 

2 1.07 X 0.038^mpiM,,^((T^u) 

Dthh — ^^2 \’ ^DM — 1/2 1/2 ’ 

gj mpi(GeV)(a^u) gj 

where rupi = \/^Mpi and g is internal degrees of freedom of DM. zdm is dehned by 
^DM = for the rjR freeze-out temperature Tf. The relevant thermally averaged 
annihilation cross section {(J-nv) including the co-annihilation processes can be found in 
[m [T5] . Since has a crucial dependence on the couplings A 3^4 given in eq. ([I]) [T5] . 

the relic abundance could change its value by varying the values of A 3^4 without 

affecting other phenomena discussed in this paper. Thus, it is not difficult to realize the 
suitable relic abundance from this source. 

The second one comes from the non-thermal origin, that is, the lepton number asym¬ 
metry left in the rj sector which is produced through the decay of A^-i. One may consider 
that this could play an important role for the DM relic abundance as in the asymmetric 
DM scenario. In fact, its contribution could be estimated as 

Ononth/i' = 2.8 X 10^1 (y^) (35) 

where is the asymmetry in the present Universe. The non-negligible contribution to 
the DM relic abundance is expected in the case AU; = 0(10“^^). However, we should note 
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that the relic abundance of rjji is fixed after the electroweak symmetry breaking. Since 
the lepton number in the rj sector is violated through the rjn-rij mass splitting caused 
by the electroweak symmetry breaking mediated by the effective coupling A 5 , the lepton 
number asymmetry in the rj sector disappears completely at this stage. Thus, this non- 
thermal component cannot contribute to the DM relic abundance in this scenario. The 
DM relic abundance is completely determined only by the thermal relics as in the same 
way discussed in the previous studies |T5]. This suggests that the leptogenesis scenario 
presented here can generate sufficient baryon number asymmetry in a consistent way with 
the generation of the neutrino masses, the DM phenomenology and others. It is notable 
that they are closely related to each other through the inflaton interaction with the SM 
Higgs scalar and 7]. 

4 Summary 

We have considered an extension of the radiative neutrino mass model with singlet scalars, 
one of which plays a role of inflaton. The original Ma model can be obtained effectively 
at low energy regions by integrating out the singlet scalars. In this model, the lepton 
number violation is prepared as the mass term of inflaton and it plays a crucial role in 
both the radiative neutrino mass generation and the generation of the lepton number 
asymmetry. The lepton number asymmetry is produced by the inflaton decay hrstly 
in the inert doublet sector. It is transferred from the inert doublet sector to the lepton 
sector through the lepton number conserving scatterings. We have examined this scenario 
numerically and showed that the sufficient baryon number asymmetry could be generated 
as long as the model parameters take suitable values. They can be consistent with the 
neutrino mass generation and the DM phenomenology. The scenario could present a new 
possibility for the leptogenesis in the framework which makes a close connection between 
the neutrino mass generation and the inflation of the Universe. 
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Appendix A 


In this Appendix, we fix the concrete form of the neutrino mass matrix to determine the 
model parameters based on the neutrino oscillation data. Since it determines the flavor 
structure of neutrino Yukawa couplings, we can fix the reaction density contained in the 
Boltzmann equations. As such a typical example, in the present analysis we use 

hgz 0; hrz — (z 1,2^, hg3 ^7-3 = /Z3, 

which could realize the tri-bimaximal neutrino mixing [9]. Although it is not realistic, it 
could give a good starting point for the purpose of this paper. In this case, three neutrino 
mass eigenvalues are given as 

= 0 , ,m^2 = 3/Z3A3, = 2 (/z^Ai + ^2^2), ( 37 ) 

where Ak is defined by 

At = Y Y A4, m,„). (38) 

a=l,2/=± ^ 

Thus, = a/A matm and = a/Atu^ should be satisfied for the normal hierarchy 
case. We use this relation to fix the values of neutrino Yukawa couplings in the present 
analysis. 


Appendix B 


In this Appendix, we give the formulas of the reaction density contributing to the Boltz¬ 
mann equations for the lepton number asymmetry. In order to give the expression for the 
reaction density of the relevant processes, we introduce dimensionless variables as 


X = 


M2' 


M 2 

xvj-r} 


^±0 


m: 


'±a 

M 2 ’ 

rj 


bu = 

/Zo 


M 2 

Tj 


where s is the squared center of mass energy. 

The reaction density for the scattering process is expressed as 

T 


'y{ab -)■ ij) = 


647r^ 


ds d{s)y/sKi 


T 


(39) 


(40) 
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where (j(s) is the reduced cross section and Ki{z) is the modihed Bessel function of the 
second kind. The lower bound of integration is dehned as Smin = max[(ma + nibY, {rrii + 
rrijf]. 

The lepton number conserving scattering processes are induced by the diagrams with 
Ni exchange which are shown in the left-hand side of Fig. 2. In order to give the ex¬ 
pression for the reaction density of these processes, we dehne the following quantities for 
convenience: 


1 X — tti 

Di{x) [x — ttiY + a^Ci 




Using these dehnitions, their reduced cross sections are expressed as 


(41) 


(Xn\X] = 


27r 




+ 


i=l 

tti 


aAx 


4a;) 


X -|- 2ai — 2 


In 


GiX + (Oj — 1)2 

X + {x'^ — 4a;) -|- 2aj — 2 


X 


(a;2 — 4a;) _|_ 2aj 


E 

i>j 


Re[{hh^)‘fj]^aiaj f 2a; -f Sa* -|- — 4 f x + (a;^ — 4a;)^/^ -|- 20 * — 2 


a; -f Oj -|- a* — 2 


Ct j 


2a; -f Oj -|- 3aj 


In 


a,- — a* 


X 


In 

a; — (a;2 — 4a;) _|_ 2aj — 2 

-f (a;^ — 4x)^/^ -I- 2aj — 2 


X — (a;2 — 4x)^/2 _|_ 2aj 


for rjr] and 


(Tb(x) = 


1 (x - 1 )^ 
27r x 2 
„2 


X 


(x-1)^ 


1 + 


X -|- a* — 2 


X 


X (x — 1)^ 


X [xa* —1 Di{x) 2Zl*(x)2 
x(x -f a* — 2) 


DAx) 


In 


xa,- — 1 




i>j 


X X (x — 1)^ 

+ + 


X 


{x-iy 


X" 


(x — 1)^ 


X IA( x) A (x) DAx)Dj{x) 

2{x + ai — 2) x-t-Oj — 2^j x(x -f a* — 2) 

xa* — 1 
x(x -f aj — 2) 


aj — a* 


Dj{x) 

2(x + aj — 2) X + ttj — 2 


a,' — a,- 


A(x) 


In 


Xttj — 1 


(42) 


(43) 


for —)■ ijsr]. 

The lepton number violating scattering processes are brought about by the diagrams 
with Fia exchange which are shown in the right-hand side of Fig. 2. In order to represent 
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their reduced cross section, we introduce the dehnition such as 


(x - h±aY + &Lc±a ’ 

2(1 - h±a) - X 






P+rv = 


Q±a — + 


2(1 - xh±a) 


1 - 


[x(a;-4)]V2 ’ 

Using these quantities, the reduced cross sections are represented as 


(Tt X = 


E "" 


dyr {x^{x — 4 )yP 


2 2 
+ 


+ 


4P_ 


Pl^-PlJ P+«-l 


In ■ 


Pl^ - 1 P 2 « - 1 

P+a + 1 


4P 


+a 


p p2 _ p2 

J^—a +0 ^—a 


In 


P-a + 1 
P-« - 1 


(cross terms between a = 1 and 2 ) 




( 44 ) 


(45) 


for ? 7?7 —)■ (pcf) and 


ay{x) 


bl 


fJ^OL 


+ 


,£^2 I Q+a - 1 Q-c. - 1 


In 


Q+a + 1 , Q-a + 1 


+ 


Q+q: Q—ol \ Q+ 

(x - 1)2 f 1 


In 


+ 


Q-a -I 
2 


4x2 ip+^(a;) P-„(x) &+«-&-« VP+«(x) D-a{x) 


X — 6+q X — 


2 x V P 


X 


-'+Q: 


X b—Qi\ f Q+a P 1 , Q—a P 1 


'^+a 


X 


In 


In 


DJx) / V Q+Q “ 1 Q-a - 1 


P (cross terms between a = 1 and 2) 


(46) 


for rjc/)'^ rf(l). Since we consider the case 3> b^^, we can neglect contributions relevant 
to bl, 


^+ 1 - 
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